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\D ! Abstract 



Modern optimisation algorithms are often metaheuristic, and they are very promising in solv- 
^) ■ ing NP-hard optimization problems. In this paper, we show how to use the recently developed 

Firefly Algorithm to solve nonlinear design problems. For the standard pressure vessel design 
optimisation, the optimal solution found by FA is far better than the best solution obtained 
43 ' previously in literature. In addition, we also propose a few new test functions with either sin- 

gularity or stochastic components but with known global optimality, and thus they can be used 
to validate new optimisation algorithms. Possible topics for further research are also discussed. 

To cite this paper as follows: Yang, X. S., (2010) 'Firefly Algorithm, Stochastic Test Functions 
and Design Optimisation', Int. J. Bio-Inspired Computation, Vol. 2, No. 2, pp. 78-84- 
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1 Introduction 

Most optimization problems in engineering are nonlinear with many constraints. Consequently, to 
. find optimal solutions to such nonlinear problems requires efficient optimisation algorithms (Deb 

1995, Baeck et al 1997, Yang 2005). In general, optimisation algorithms can be classified into two 
main categories: deterministic and stochastic. Deterministic algorithms such as hill-climbing will 
produce the same set of solutions if the iterations start with the same initial guess. On the other 
hand, stochastic algorithms often produce different solutions even with the same initial starting 
point. However, the final results, though slightly different, will usually converge to the same optimal 
solutions within a given accuracy. 

Deterministic algorithms are almost all local search algorithms, and they are quite efficient in 
finding local optima. However, there is a risk for the algorithms to be trapped at local optima, while 
the global optima are out of reach. A common practice is to introduce some stochastic component 
to an algorithm so that it becomes possible to jump out of such locality. In this case, algorithms 
become stochastic. 

Stochastic algorithms often have a deterministic component and a random component. The 
stochastic component can take many forms such as simple randomization by randomly sampling the 
search space or by random walks. Most stochastic algorithms can be considered as metaheuristic, 
and good examples are genetic algorithms (GA) (Holland 1976, Goldberg 1989) and particle swarm 
optimisation (PSO) (Kennedy and Eberhart 1995, Kennedy et al 2001). Many modern metaheuristic 
algorithms were developed based on the swarm intelligence in nature (Kennedy and Eberhart 1995, 
Dorigo and Stiitzle 2004). New modern metaheuristic algorithms are being developed and begin to 
show their power and efficiency. For example, the Firefly Algorithm developed by the author shows 
its superiority over some traditional algorithms (Yang 2009, Yang 2009, Lukasik and Zak 2009). 

The paper is organized as follows: we will first briefly outline the main idea of the Firefly 
Algorithm in Section 2, and we then describe a few new test functions with singularity and/or 
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randomness in Section 3. In Section 4, we will use FA to find the optimal solution of a pressure 
vessel design problem. Finally, we will discuss the topics for further studies. 

2 Firefly Algorithm and its Implementation 

2.1 Firefly Algorithm 

The Firefly Algorithm was developed by the author (Yang 2008, Yang 2009), and it was based on 
the idealized behaviour of the flashing characteristics of fireflies. For simplicity, we can idealize these 
flashing characteristics as the following three rules 

• all fireflies are unisex so that one firefly is attracted to other fireflies regardless of their sex; 

• Attractiveness is proportional to their brightness, thus for any two flashing fireflies, the less 
brighter one will move towards the brighter one. The attractiveness is proportional to the 
brightness and they both decrease as their distance increases. If no one is brighter than a 
particular firefly, it moves randomly; 

• The brightness or light intensity of a firefly is affected or determined by the landscape of the 
objective function to be optimised. 

For a maximization problem, the brightness can simply be proportional to the objective function. 
Other forms of brightness can be defined in a similar way to the fitness function in genetic algorithms 
or the bacterial foraging algorithm (BFA) (Gazi and Passino 2004). 

In the FA, there are two important issues: the variation of light intensity and formulation of the 
attractiveness. For simplicity, we can always assume that the attractiveness of a firefly is determined 
by its brightness or light intensity which in turn is associated with the encoded objective function. 
In the simplest case for maximum optimization problems, the brightness I of a firefly at a particular 
location x can be chosen as / (x) oc /(x). However, the attractiveness ft is relative, it should be seen 
in the eyes of the beholder or judged by the other fireflies. Thus, it should vary with the distance 
rij between firefly i and firefly j. As light intensity decreases with the distance from its source, and 
light is also absorbed in the media, so we should allow the attractiveness to vary with the degree of 
absorption. 

In the simplest form, the light intensity I(r) varies with the distance r monotonically and expo- 
nentially. That is 

J = Ioe-t^ (1) 

where Io is the original light intensity and 7 is the light absorption coefficient. As a firefly's at- 
tractiveness is proportional to the light intensity seen by adjacent fireflies, we can now define the 
attractiveness (3 of a firefly by 

P = O e~^\ (2) 

where /?o is the attractiveness at r = 0. It is worth pointing out that the exponent can be 
replaced by other functions such as jr m when m > 0. Schematically, the Firefly Algorithm (FA) 
can be summarised as the pseudo code. 
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Firefly Algorithm 

Objective function /(x), x= (x\, Xd) T 
Initialize a population of fireflies Xj (i = 1, 2, n) 
Define light absorption coefficient 7 
while (t <MaxGeneration) 
for i = 1 : n all n fireflies 
for j = 1 : i all n fireflies 

Light intensity 7j at x, is determined by /(xj) 
if (/,->/<) 

Move firefly i towards j in all d dimensions 
end if 

Attractiveness varies with distance r via exp[— -yr] 
Evaluate new solutions and update light intensity 
end for j 
end for i 

Rank the fireflies and find the current best 
end while 

Postprocess results and visualization 

The distance between any two fireflies i and j at Xj and x.,- can be the Cartesian distance 
Vij = ||xj — 1 1 2 or the ^ 2 -norm. For other applications such as scheduling, the distance can be time 
delay or any suitable forms, not necessarily the Cartesian distance. 

The movement of a firefly i is attracted to another more attractive (brighter) firefly j is deter- 
mined by 

Xj = Xj + /3 e~ 7r ^ (xj - Xj) + aei, (3) 

where the second term is due to the attraction, while the third term is randomization with the vector 
of random variables being drawn from a Gaussian distribution. 

For most cases in our implementation, we can take ft = 1, a e [0, 1], and 7 = 1. In addition, 
if the scales vary significantly in different dimensions such as — 10 5 to 10 5 in one dimension while, 
say, — 10~ 3 to 10 3 along others, it is a good idea to replace a by aSk where the scaling parameters 
Sk(k = l,...,d) in the d dimensions should be determined by the actual scales of the problem of 
interest. 

In essence, the parameter 7 characterizes the variation of the attractiveness, and its value is 
crucially important in determining the speed of the convergence and how the FA algorithm behaves. 
In theory, 7 € [0, 00), but in practice, 7 = 0(1) is determined by the characteristic/mean length Sk 
of the system to be optimized. In one extreme when 7 — >■ 0, the attractiveness is constant j3 — 
This is equivalent to saying that the light intensity does not decrease in an idealized sky. Thus, a 
flashing firefly can be seen anywhere in the domain. Thus, a single (usually global) optimum can 
easily be reached. This corresponds to a special case of particle swarm optimization (PSO). In fact, 
if the inner loop for j is removed and Ij is replaced by the current global best g*, FA essentially 
becomes the standard PSO, and, subsequently, the efficiency of this special case is the same as that 
of PSO. On the other hand, if 7 — > 00, we have j3{r) — > S(r), which is a Dirac 5-function. This means 
that the attractiveness is almost zero in the sight of other fireflies, or the fireflies are short-sighted. 
This is equivalent to the case where the fireflies fly in a very foggy region randomly. No other 
fireflies can be seen, and each firefly roams in a completely random way. Therefore, this corresponds 
to the completely random search method. So 7 partly controls how the algorithm behaves. It is also 
possible to adjust 7 so that multiple optima can be found at the same during iterations. 

2.2 Numerical Examples 

From the pseudo code, it is relatively straightforward to implement the Firefly Algorithm using a 
popular programming language such as Matlab. We have tested it against more than a dozen test 
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Figure 1: Four global maxima at (±1/2, ±1/2). 



functions such as the Ackley function 



/(x) 



) = -20cx P [-i 

^ i=i 



exp[- cos(27ra; i )] + 20 + e, 



(4) 



which has a unique global minimum /* = at (0,0,..., 0). From a simple parameter studies, we 
concluded that, in our simulations, we can use the following values of parameters a. = 0.2, 7 = 1, 
and po = 1. As an example, we now use the FA to find the global maxima of the following function 



with the domain —10 < Xi < 10 for all (i = 1, 2, ...,d) where d is the number of dimensions. This 
function has multiple global optima. In the case of d = 2, we have 4 equal maxima /* = ~ 



0.6065 at (1/2,1/2), (1/2,-1/2), (-1/2,1/2) and (-1/2,-1/2) and a unique global minimum at 
(0,0). 



The four peaks are shown in Fig. 1, and these global maxima can be found using the implemented 
Firefly Algorithms after about 500 function evaluations. This corresponds to 25 fireflies evolving 
for 20 generations or iterations. The initial locations of 25 fireflies are shown Fig. 2 and their 
final locations after 20 iterations are shown in Fig. 3. We can see that the Firefly Algorithm is 
very efficient. Recently studies also confirmed its promising power in solving nonlinear constrained 
optimization tasks (Yang 2009, Lukasik and Zak 2009). 



The literature about test functions is vast, often with different collections of test functions for 
validating new optimisation algorithms. Test functions such as Rosenbrock's banana function and 
Ackley's function mentioned earlier are well-known in the optimisation literature. Almost all these 
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3 New Test Functions 
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Figure 2: Initial locations of 25 fireflies. 




Figure 3: Final locations after 20 iterations. 
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test functions are deterministic and smooth. In the rest of this paper, we first propose a few new 
test functions which have some singularity and/or stochastic components. Some of the formulated 
functions have stochastic components but their global optima are deterministic. Then, we will use 
the Firefly Algorithm to find the optimal solutions of some of these new functions. 
The first test function we have designed is a multimodal nonlinear function 

/(x)=[ e -Et 1 (^)-_ 2e -Et 1 (— ) 2 

d 

•J^Jcos 2 Xi, m = 5, (6) 

i=l 

which looks like a standing- wave function with a defect (see Fig. 4). It has many local minima 
and the unique global minimum /* = —1 at x* = (it, tt, ir) for j3 = 15 within the domain 
—20 < Xi < 20 for i = 1,2, d. By using the Firefly Algorithm with 20 fireflies, it is easy to find 
the global minimum in just about 15 iterations. The results are shown in Fig. 5 and Fig. 6. 

As most test functions are smooth, the next function we have formulated is also multimodal but 
it has a singularity 



/(x) = (^|x,|).cxp[-^sin(x 2 )], (7) 



which has a unique global minimum /* = at x* = (0, 0, 0) in the domain — 2ir < Xj < 2ir where 
i = 1,2, ...,d. At a first look, this function has some similarity with function (5) discussed earlier. 
However, this function is not smooth, and its derivatives are not well defined at the optimum 
(0,0,..., 0). The landscape of this forest-like function is shown in Fig. 7 and its 2D contour is 
displayed in Fig. 8. 

Almost all existing test functions are deterministic. Now let us design a test function with 
stochastic components 

f(x,y) = -5e-W x -^ 2+ (y-^ 

-EX>.? e ~ a[(M)2+( ^A (8) 

j=l i=l 

where a, (3 > are scaling parameters, which can often be taken as a = /3 = 1. Here the random 
variables (i, j — 1, K) obey a uniform distribution ~ Unif[0,l]. The domain is < x, y < K 
and K = 10. This function has K 2 local valleys at grid locations and the fixed global minimum at 
x» = (7r,7r). It is worth pointing that the minimum / m j n is random, rather than a fixed value; it 
may vary from —{K 2 + 5) to —5, depending a and (3 as well as the random numbers drawn. 

For stochastic test functions, most deterministic algorithms such as hill-climbing would simply 
fail due to the fact that the landscape is constantly changing. However, metaheuristic algorithms 
could still be robust in dealing with such functions. The landscape of a realization of this stochastic 
function is shown in Fig. 9. 

Using the Firefly Algorithm, we can find the global minimum in about 15 iterations for n = 20 
fireflies. That is, the total number of function evaluations is just 300. This is indeed very efficient 
and robust. The initial locations of the fireflies are shown in Fig. 10 and the final results are shown 
Fig. 11. 

Furthermore, we can also design a relative generic stochastic function which is both stochastic 
and non-smooth 

d 

/(x) = 22 £i \Xi , -5 < x t < 5, (9) 
i=i 

where (i = l,2,...,d) are random variables which are uniformly distributed in [0,1]. That is, 
£i ^Unif[0, 1]. This function has the unique minimum /* = at x„ = (0,0,..., 0) which is also 
singular. 

We found that for most problems n = 15 to 50 would be sufficient. For tougher problems, larger 
n can be used, though excessively large n should not be used unless there is no better alternative, 
as it is more computationally extensive for large n. 
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Figure 4: The standing wave function for two independent variables with the global minimum 
/* = -1 at (7T,7r). 
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Figure 6: Final locations after 15 iterations. 
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Figure 7: The landscape of function (7). 
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Figure 8: Contour of function (7). 









Figure 9: The 2D Stochastic function for K = 10 with a unique global minimum at (7r,7r), though 
the value of this global minimum is somewhat stochastic. 
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Figure 10: The initial locations of 20 fireflies. 




Figure 11: The final locations of 20 fireflies after 15 iterations, converging into (n, n). 
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4 Engineering Applications 



The simple bounds are 
and 



Now we can apply the Firefly Algorithm to carry out various design optimisation tasks. In principle, 
any optimization problems that can be solved by genetic algorithms and particle swarm optimisation 
can also be solved by the Firefly Algorithm. For simplicity to demonstrate its effectiveness in real- 
world optimisation, we use the FA to find the optimal solution of the standard but quite difficult 
pressure design optimisation problem. 

Pressure vessels are literally everywhere such as champagne bottles and gas tanks. For a given 
volume and working pressure, the basic aim of designing a cylindrical vessel is to minimize the total 
cost. Typically, the design variables are the thickness d\ of the head, the thickness d 2 of the body, 
the inner radius r, and the length L of the cylindrical section (Coello 2000, Cagnina et al 2008). 
This is a well-known test problem for optimization and it can be written as 

minimize /(x) = 0.6224dir£ + 1.7781d 2 r 2 

+3.1661d?L+19.84d?r, (10) 

subject to the following constraints 

gi (x) = -di + 0.0193r < 

02 (x) = -da + 0.00954r < , , 

g 3 (x) = -irr 2 L- ^-r 3 + 1296000 < [ ' 

04 (x) = L - 240 < 0. 

0.0625 < di, d 2 < 99 x 0.0625, (12) 

10.0 <r, L < 200.0. (13) 

Recently, Cagnina et al (2008) used an efficient particle swarm optimiser to solve this problem 
and they found the best solution 

/* w 6059.714, (14) 

at 

x* w (0.8125, 0.4375, 42.0984, 176.6366). (15) 

This means the lowest price is about $6059.71. 

Using the Firefly Algorithm, we have found an even better solution with 40 fireflies after 20 
iterations, and we have obtained 

x, w (0.7782, 0.3846, 40.3196, 200.0000) T , (16) 

with 

/ min « 5885.33, (17) 

which is significantly lower or cheaper than the solution /* w 6059.714 obtained by Cagnina et al 
(2008). 

This clearly shows how efficient and effective the Firefly Algorithm could be. Obviously, further 
applications arc highly needed to see how it may behave for solving various tough engineering 
optimistion problems. 

5 Conclusions 

We have successfully used the Firefly Algorithm to carry out nonlinear design optimisation. We first 
validated the algorithms using some standard test functions. After designing some new test functions 
with singularity and stochastic components, we then used the FA to solve these unconstrained 
stochastic functions. We also applied it to find a better global solution to the pressure vessel design 
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optimisation. The optimisation results imply that the Firefly Algorithm is potentially more powerful 
than other existing algorithms such as particle swarm optimisation. 

The convergence analysis of metaheuristic algorithms still requires some theoretical framework. 
At the moment, it still lacks of a general framework for such analysis. Fortunately, various studies 
started to propose a feasible measure for comparing algorithm performance. For example, Shilane et 
al (2008) suggested a framework for evaluating statistical performance of evolutionary algorithms. 
Obviously, more comparison studies are highly needed so as to identify the strength and weakness 
of current metaheuristic algorithms. Ultimately, even better optimisation algorithms may emerge. 
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